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Abstract— This paper deals with design optimization of 
asynchronous machine (3 HP, 4-POLE, 3-PHASE, 50 Hz) 
considering three objective functions torque (T), cost (C) and 
efficiency (n); using genetic algorithm (GA) / improved genetic 
algorithm (IGA) & non-linear programming. 
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I. INTRODUCTION 


Induction motors are being widely used for almost all 
practical purposes in industrial & domestic world, giving rise 
to the consumption of maximum share of the total electrical 
power generated. So, their efficiency, torque produced and cost 
becomes a major concern to be optimized. But the optimal 
design of I.M. is a problem due to the following reasons: 


1. Induction motor involves so many variables, which non- 
linearly affects the performance of machine. 

2. Due to involvement of various conflicting parameters, 
objectives and their non-linear behavior during design and 
operation of the machine, it becomes quite necessary to 
balance between many variables, sacrificing one for the 
other to obtain the best possible performance and 
optimized design depending on end user’s application. 

The effect of motor performance on surroundings and 
environment is of great concern also, e.g., the level of noise and 
acoustic comfort of train passengers and nearby residents has to 
be ensured. The motor starts making noise due to Maxwell's air 
gap forces, which causes in creating stator vibrations in audible 
range. 


The problem of induction motor optimal design has received 
much attention since the beginning of computer sciences 
[4,13,16]. Many solutions and algorithms were designed based 
on following points. 


Actually many conventions are made on the basis of 
practical experiences with manipulations of Induction motor 
design, e.g., a small air gap improves the efficiency, but 
increases the magnetic sound level and decreases the overload 
capacity of the motor; increasing the stator length of yoke to 
diameter ratio generally lowers magnetic vibrations [14]; but it 
increases the material cost; and with a fixed motor size, it 
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decreases the available out-put torque with increased rotor 
temperature. 


In literature so many single as well as multiple objective 
approaches for the electrical machine design has been proposed 
[13,17]; using various methods and techniques. These multiple 
objective approaches may be dealt with several NLP methods 
[8,9,10,11], as well as genetic algorithm [5,6,7,12,18] and 
improved genetic algorithm [15]. 


II. PROBLEM DEFINITION & DESIGN APPROACH 


The standard single phase equivalent circuit model of a 3- 
phase induction motor on per phase basis is shown in figure-1. 
The model offers reasonably good prediction accuracy with 
modest computational efforts, despite its shortcomings. This 
model is basically a per phase representation of a balanced poly- 
phase induction motor in the frequency domain, comprising six 
model parameters. The six parameters have their usual meaning 


[5]. 


III, AN OVERVIEW OF OPTIMIZATION BY GENETIC ALGORITHM 


In the most general sense, GA-based optimization is a 
stochastic search method that involves the random generation of 
potential design solutions and then systematically evaluates and 
refines the solutions until a stopping criterion is met. There are 
three fundamental operators involved in the search process of a 
genetic algorithm: selection, crossover, and mutation. The 
genetic algorithm implementation steps are shown as follows: 


1. Define 
(Initializing). 

2. Generate first population in a random manner from 
search space. 


parameters and objective _ functions 


Ro'(1/s-1) 


Fig. 1 Equivalent circuit model of an induction motor 


50 


3. Evaluate population by using objective functions and 
@srange in order of merit. 

4. Test convergence. If satisfied then stop else continue. 

5. Start reproduction process (Selection, Crossover, 
Mutation & elitism). 

6. Form new generation of offspring and treat as new 
population. To continue the optimization, return to step 3. 

To apply the GA approach, objective functions and 
constraints have to be defined to evaluate how good the motor 
design is obtained. 


IV. IMPLEMENTATION OF THE OPTIMAL DESIGN PROCEDURE 


The formulation of the objective function [5] is as per the 
following process: 


F(x) =} FOr) - Plg; 71, 
0, If F(x)- Plg,(x),r]s0 


If F(x) — P[g;(x),r] > 0 
(1) 


Where, F(x) is the objective function, motor material 
cost/efficiency/torque, r is the penalty coefficient related to the 
value of objective function and x is design variable set. The 
penalty functions, P [g,(x), r], are expressed with respect to the 
type of inequality used. 


By means of exterior penalty function, constrained problems 
are converted to unconstrained problems by removing 
constraints. According to constraints, penalty function is 
defined as following [5]: 


P(g,(x),r) =r [max (0,9,)], § = 1,2,..m 


r [min (0,g;)]7, <i 


(2) 


When the constraint inequality is satisfied, the penalty 
function becomes inactive. The objective function emphasizes 
the larger constraint violations and the optimization search 
tries to reduce these violations to zero, in the feasible region. 
This would result in pushing the search into the feasible 
design region. All the constraints are satisfied within this 
region and the optimization approach attempts to move the 
design into its best optimum solution [5]. 


j=m,. 


However, order of magnitudes of various constraints is 
much different from one another in case of electric motors. 
Therefore, constraint functions need to be normalized with 
respect to the specified objective function to have meaningful 
convergence criteria, It is to be ensured that constraints with 
higher values do not dominate over others. The normalized 
constraint functions in the penalty function are developed as 
shown in the following [5]. 


Dj norm(X) = Birer — 2))/bp f=1,2,...0 (3) 


Where bjrer is the calculated value from the current 
generation / evaluation whereas bj is the expert defined 
constraint as shown in Table 1. The main purpose for defining 
the constraint b; is to have the final design for practically be 
feasible and acceptable. In general, the constraints are decided 
upon with great care taking into consideration the availability 
of materials, customer’s requirements and manufacturing 
standards. Table 1 is also referred to as the motor 
specifications and their constraint values. Constraint values of 
variables can be expressed by following inequality [5] given 
below. 


9G) =s—b,< 0, ge(x) = pe 
92 (x) = By-b; <0, 9;(x) = Cos? — b, >0 
T, 
93) = Byy-b3 <0, ax) => 20 
— bs 
9.0%) = By, — 6, <0, 99(X) = =_59 
T—b, 


95x) = Fy — bs <0, G10(X) = 1 — by 
(4) 

In eqn.4, there are two different conditions of inequality 

constraint which are explained in following sections. 


TABLE I INEQUALITY TABLE FOR MOTOR PARAMETERS 


SR.NO, NAME OF PARAMETER INEQUALITY (bj) 

I Slip. s = (bi = 0.05) 

2 Stator yoke flux density, By < (b= 0.6) 

3 Rotor yoke flux density, B, < (b; = 0.6) 

4 Stator teeth flux density, Be < (bs = 1.7) 

5 Stator slot filling factor, F; < (bs = 0.90) 
Starting current to rated 

6 current ratio, Ijun / I S (bs =7) 

[7 Power factor, cos® 2 (by = 0.8) 

Pull-out torque to rated 

8 torque, T, / T 2 (bs= 1.5) 
Starting torque to rated 

9 torque ratio T/T = (by = 0.6) 

10 Efficiency (n) > (bio = 0.85) _| 


A. First Condition 


It is not permitted that some constraints [as g,(x), j=7.... 10] 
are below the expert defined level. For example, high value of 
power factor is desired for good performance in induction 
motors. If the expert-defined constraint for power factor as 
shown in Table | was 0.8, then anything less than that would 
be a violation. 


B. Second Condition 


It is not permitted that some constraints [as g(x), j=1,..6] 
are above the expert defined level. For example, stator yoke 
flux density may not exceed certain values on account of iron 
losses. If the expert-defined constraints for stator yoke flux 
density as shown in Table 1 was 0.6, then anything more than 
that would be a violation. 
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Fig. 2 Flow chart for 3-phase induction motor design using NLP 
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Fig. 3 flow chart for optimization usina GA/ IGA 

The software developed for the design optimization of the 
induction motor was prepared using “JAVA”, which can 
analyze, optimize & evaluate design parameters and 
performance of motor. Parameters of the motor or materials 
used can be easily modified to investigate their effect on 
performance. Selection and optimization type (efficiency, 
torque, cost etc.) can be made depending on user. 


The GA optimization algorithm was based on a roulette 
wheel selection, single point crossover, bit mutation and 
elitism. The flow chart for design optimization process using 
NLP is given in figure 2 and the flow chart for design 
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Fig. 4 Flow-chart for induction 


design using GA/IGA 
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optimization process using GA is given in Fig. 3 & 4. Each 
block consists of a number of subroutines. 


V. MODEL FOR 3-PHASE INDUCTION MOTOR DESIGN V. 
OPTIMIZATION 


The highly nonlinear constrained ~—_—multivariable 
optimization problem is very difficult to be solved by the 
conventional methods. The optimization problem for 3-phase 
induction motor design can be formulated as [15]: 


Min f(x) 
glx) 20, b= 1:2;3,.0.5n 
A(x)20 é = 1,2,3,....,n 


(5) 


Where f(x) is the objective function of optimization, gi(x) 
and h(x) are constraining functions, and x is the design 
variable set. 


A. Objective Functions 


In order to reduce the active (main) material cost and 
improve the efficiency of 3-phase induction motors, three 
different objective functions of optimization are defined 
separately as: [15] 


f, (x) = Ce, p uWse = Coy pu Weu sr Cy pu) ro 


A£@)=n A(x) =T (6) 
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Where Crepus Ccupus Caipu are the unit prices of magnetic 


material, copper wire and aluminum conductive _ bars 
respectively, Wre, Wcu, and Wa; are the masses of the 
corresponding materials used in the construction of motor, 1 is 
the motor efficiency and T is full load torque produced. 


_ Pout.put 
| aan (7) 
Pout.puttLosses 
Where, Losses = Stator iron loss + Stator copper loss + Rotor 
copper loss + mechanical loss, 
- Losses = y 3V, 1, cos0- Losses 
(8) 


& P, 


output 


=P, 


input 


B. Design Optimization Variables 


There is a lot of design variables involved in the design of 
a 3-phase induction motor and it is required to choose 
critically the suitable parameters as the optimization design 
variables as per the following carefully [15]: 


x , et ep aided: 


[L,D, lg, Z,Hs,We. ly, Dy, ATM, Bay: Seu Sal 


[L, Dy lg Zp Hy My, bpp Bp ATM, Boy 


(9) 


Where L stator core length, D stator core inner diameter, |, 
length of air gap, Z number of stator conductors, H, height of 
stator slot, W, width of stator slot, |, height of rotor slot, b, 
width of rotor slot, ATm ampere conductors per meter, Bay 
average flux density, 5c, current density in stator winding 
conductor, 5,; current density in rotor bars. 


C. Constraining Functions 


The Genetic Algorithm usually adopts the sequential 
unconstrained minimization technique (SUMT) to solve the 
problems defined by previous article [(V) B]. Experiences and 
simulations have shown that if the penalty factor is selected 
too large, the CGA may cause premature convergence. On the 
other hand, if the penalty factor is selected too small, the CGA 
may yield a highly computational burden. To overcome it, a 
new annealing penalty function, which adopts self-adaptive 
annealing penalty factors, is developed for Improved genetic 
algorithm (IGA) [15]: 


(10) 
f(x,o) = f(x) + p(x, o) 
p(x,a) = o[Zi imin(O, gi(x)| + Xi thi(x)I] 


where, 


o is the penalty factor (O=1/T, Ti+1=Ta, , T is the value of 
the parameter under control, a is the proportionality constant 
within the range of 0-1), p(x, 6) the penalty function, and f(x, 
o) the generalized objective function defined by the IGA. 


To reflect the changing feature of the nature and imitate the 
parameter under control becoming more rigorous with the 
development of the generations, the fitness value of a string is 
defined as 


fit= f,,—f (Xx) (11) 


Where, f,, is a maximum objective function in a generation. 
To avoid prematurity, another technique known as the multi- 
turn evolution strategy was used. 


D. Application 
Optimization 


of IGA to 3-Phase I.M. Design 

The IGA is applied for design optimization of the 3-phase 
induction motor. In the process of optimization, the objective 
as well as constraining functions in previous articles [(IV), (V) 
C] are implemented as the following [15]: 


g3() = 


gs) 


g74x) = Ip, 2 0, 520 7 bp 2 


ggix) =ATm- ATm' = 0.92500) = Boy Bay 


9:1) = Sq. 8h, 20, 


(12) 


Symbols have their meanings as given in eqn 9. 


VI. OBJECTIVE FORMATIONS 


The major performance parameters, which affects the 
machine performance typically, can be collected as: 1) Cost of 
machine, 2) Efficiency & Torque (FL) produced, 3) Break 
down torque, 4) Starting torque, 5) Power factor of the 
machine, 6) Full load slip of machine, 7) Starting to rated 
current ratio, 8) Pullout to rated torque ratio (over load 
capacity of the motor), 9) Magnetic vibrations & noise. 


A. Cost of Machine 
It involves following parameters: 


COST OF IRON: 


Cre = L = OD* wk, * Pee * Crepu (13) 
COST OF COPPER: 

Cus = (14) 
COST OF ALUMINUM: 

Caiy = Ly * Ap * 52% Pai > Caipn (15) 


COST OF INSULATION: 


For a given voltage, Cost of Insulation material is proportional 
to Temperature rise of Machine & hence is proportional to 
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total losses of the machine, and will be considered in 
efficiency maximization. 


Total cost of the machine can be given by: 


(16) 


Ce = Cre> Cons + Coir + Cinsuiation + Cy 


Where C, is punching cost and can be taken as 20% of total 
cost. 


B. Efficiency and Torque Produced at Full Load 


The function for minimizing the losses can be used as 
function for maximizing the efficiency. 


We know that, 


Ef ficiency(n) a motte (17) 
And Losses in the motor can be given by: 
Total Losses = W, + We, + F &w loss 
(18) 


For maximum efficiency of the motor at full load, iron losses 
can be taken as equal to Poy stator; friction and windage loss can 
be taken as 10% of full load rotor copper loss (Pex rotor). 


So, now total losses can be given by: 


Total Losses = 2Pex stator + 1. (19) 

Total Losses = 2 » 313,Revag + 

1.112 (RS; +Rring) 
(20) 
C. Breakdown (Pull-out) Torque 
1 

Tgp @& | ---+——_ 

BD rotor reactance(Xy) Ql) 


D. Starting Torque 


We can obtain better starting torque with R / X ratio near 
to one at standstill, i.e. 


If R,/%, =1: 


Better starting torque can also be obtained by deep bar or 
double cage rotor. 


Tye a R,: T,; = Tsp 


E. Starting Torque 


We know very well that power factor of the motor is better 
with smaller air gap length (1p). _ 


pra (7) 
: (22) 
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F. Full Load Slip 


= Pou 
We know that slip, Pg (23) 
or, s a Copper Losses (24) 


G. — Starting to Rated Current Ratio 


Starting Current is low with high rotor resistance, so it can 
be understood that if R / X ratio is near to one for better 
starting torque, the starting current will be low and ratio of 
starting current to rated current will be better for such a motor, 
but at the cost of efficiency. 


H.  Pull-out to Rated Torque Ratio 
Pull-out to rated torque ratio 


= (Teuiteus) t ‘ 
( Traore } ® am mae 


(25) 


I, Magnetic Vibrations and Noise 

We know that noise caused due to magnetic vibrations is a big 
problem, which is due to the zig-zag leakage reactance, which 
is inversely proportional to the air-gap length of the motor, i.e. 
Noise a (Zig — Zag leakage reactance) (26) 
But on the contrary larger air gap give rise to the magnetizing 
reactance which causes lower power factor and hence inferior 
performance of the motor. So, one has to make sacrifice 
between the noise and power factor of the motor. 


VU. PERFORMANCE INDEX 


Now we can write performance index as: 


(27) 


Where, kj, k»....ks are weight constants given to different 
functions, given in table 2. 


TABLE II WEIGHTS FOR PERFORMANCE INDEX 


NAME OF FUNCTION WEIGHT GIVEN TO 
THE FUNCTION 


|__1__| Cost of the motor for 


2 Total losses 
{3 Starting torque 
4 F (Length of air gap, lg) 


[3s] Palo ote org a [or 


Dee os 
Function 
ad 
Index 


These weights are decided for a motor running for a small 
Atta-Chakki under following conditions: 


1. Continuous running for more than 20 hrs per day. 


ie) 


we 


The motor should be economic but due to its continuous 


running and smaller size (lower capital cost) energy 
efficiency and performance is more important. 


The motor is to be used in a silent environment (say 
residential area, or for house hold use). So, we have to 
consider weight for noise function also, but it is very less 
due to the fact that it conflicts with the power factor and 
hence the efficiency of the motor. One has to sacrifice 
between the noise and the efficiency and performance of 
the motor. 


TABLE II] PARAMETERS USED IN GA/ IGA 


SR.NO. | PARAMETER VALUE / CRITERIA 

1 Population Size 400 

2 Cross over rate 0.8 

3 Mutation rate | 0.01 

4 Elite count 2 
|Performance Index| < 0.15 
&E& 

5 Stopping criteria {1 / {| Performance Index_1 - 
Performance Index_0 |}] = 
null 

6 | Number of generations 10,000 (Max) 


Single point crossover with 


7 Cinsn over niettiod fixed cross over rate (0.8) 
A Bit string mutation with fixed 
8 Mutation method mutation rate (0.01) 
9 Selection coefficient Fitness 
10 Fitness method Genetic load 
VIII. RESULTS 


Results are summarized and analyzed in Table 4. These 


results are also depicted in Figs. 5, 6 & 7 given at the end. 


IX. CONCLUSIONS 


Result analysis of design, optimization and validation (by 


MATLAB/SIMULINK MODEL) gives 


the following 


conclusions: 


1 


1] 
2] 
B] 
[4] 


(5] 


The output of the GA program and its validation using 
SIMULINK / MATLAB shows that: 


The cost of motor decreases with increase in efficiency. 


Starting current of motor increases considerably but with 
better starting torque, which is required in the assumed 
case, as on load starting is required. Breakdown torque is 
also improved considerably and hence overload capacity 
also increases. 


The optimized motor design using GA/IGA reaches full 
speed faster having lesser transient time, lower slip & 
higher efficiency at full load and lower cost, so the 
temperature rise of the windings and hence insulation 
during starting as well as running will be lower and the 
life of the motor will be longer due to longer life of the 
insulation material. 
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